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The thermal Wightman functions for free, massless parti- 
cles of spin 0, 1/2, 1, 3/2, and 2 are computed directly in coor- 
dinate space by solving the appropriate differential equation 
and imposing the Kubo-Martin-Schwinger condition. The so- 
lutions are valid for real, imaginary, or complex time. The 
Wightman functions for spin 1 gauge bosons and for spin 2 
gravitons are directly related to the fundamental functions 
for spin 0. The Wightman functions for spin 3/2 gravitinos 
is directly related to that for spin 1/2 fermions. Calculations 
for spin 1, 3/2, and 2 are done in covariant gauges. In the 
deep space-like region the Wightman functions for bosons fall 
like T/r whereas those for the fermions fall exponentially. In 
the deep time-like region all the Wightman functions fall ex- 
ponentially. 
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I. INTRODUCTION 

In both zero-temperature and finite-temperature field 
theory it is customary to perform calculations in momen- 
tum space. There are exceptions to this pattern however. 
Lattice gauge theory computations are done in coordi- 
nate space. The short-distance operator product expan- 
sion is formulated in coordinate space. It is even possible 
to carry out ultraviolet regularization and renormaliza- 
tion directly in coordinate space m . 

The purpose of this paper is to deduce the free Wight- 
man functions at finite temperature for various massless 
particles directly in coordinate space. Thermal averages 
are performed with respect to the equilibrium density 
operator 
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where H is the appropriate Hamiltonian and [3 — l/T. 
Let (p (x) denote a quantum field of any spin. For spin- 
zero fields the index A only distinguishes species. For 
spin 1/2 fields A denotes a spinor index; for spin 1, a 
vector index; for spin 3/2 a spinor and vector index; for 
spin 2, a pair of vector indices. The Wightman functions 
at finite temperature are 



G^^(x) = -zTr(e0^(x)0^(O)) 
G^^(x) = -*Tr(g0^(O)0^(a;))(-l)2''. 



(1.2) 



Knowing the thermal Wightman function allows direct 
construction of the various thermal propagators in real 
or imaginary time M-m- In particular the time-ordered 
Green function is 



The canonical method for obtaining G> {x) for free 
fields would require the following steps: (1) Solve the free 
field equations and express the field operator in terms 
of plane-wave solutions weighted by creation and anni- 
hilation operators. (2) Impose the equal-time canonical 
commutation relations. (3) Express the Hamiltonian in 
terms of creation and annihilation operators. (4) Com- 
pute the trace in Fock space over the density operator. 
(5) Perform the integral over momentum states so as to 
obtain G> (x) in coordinate space. For higher spins this 
becomes tedious. 

The first four steps can be avoided if the T = prop- 
agator is already known, because the spectral function 
(K) gives directly the Wightman function in momen- 
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tum space (^-^: 
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However it is still necessary to Fourier transform from 
momentum space to coordinate space and obtain G>(a;). 
This was the procedure followed in [p|, which computed 
the coordinate-space Wightman functions for massless 
vector bosons in various gauges (Feynman, general co- 
variant, and Coulomb). However performing the Fourier 
transforms so as to ensure the correct analyticity prop- 
erties in complex time is difficult. For gauge bosons the 
complete answer satisfying all the analyticity properties 
was only obtained in the Feynman gauge. It would be la- 
borious to pursue the Fourier transform method for spin 
1/2 fermions, for gravitinos, or for gravitons. 

It turns out to be simpler to deduce the thermal Wight- 
man functions for free fields by working directly in coor- 
dinate space. Since G>(a;) solves the free field equation 
all that is necessary is that the solution have the cor- 
rect zero-temperature limit and satisfy the Kubo-Martin- 
Schwinger periodicity relation 0-5,01 under t —^ t — i[3. It 
has not been generally recognized that the KMS relation 
is not only a necessary condition that thermal Wightman 
functions must satisfy but is also sufficient condition to 
determine them directly. 

The paper is organized according to spin. Sec II deals 
with bosons. The thermal Wightman functions for spin- 
l ess b osons is given in Eq. (2JX]|) ^ for gauge boso ns in Eq. 
( ^.171) , and for gravitons in Eq. (|2.22| ) and (|2.28D . Sec HI 
deals with fermions. The thermal W ight man functions 
for spin 1/2 fermions is given by Eq. (3.3) and (3.S); for 
spin 3/2 gravitinos in Eq. (3.12) and ( |3.14 ). In the co- 
variant gauges considered here, the higher spin functions 



are all expressible in terms of the basic spin and 1/2 
functions. Sec III discusses the asymptotic behavior. 

Throughout the paper the fundamental functions de- 
pend only radial distance r and time t and it is convenient 
to use variables u and v defined as follows: 



r + t 



t. 



(1.5) 



To avoid confusion with the metric g^i, in curved space, 
the Minkowski metric will everywhere be denoted by rj^^. 



II. BOSONS OF SPIN 0,1,2 

A. Spinless Bosons 

For a spinless boson field (i>{x) the basic thermal Wight- 
man functions are 



i?>(x) = -*Tr(ec/.(a:)0(O)) 
i?<(x) = -iTr(e0(O)0(a;)). 



(2.1) 



The emphasis of the subsequent development will be to 
avoid expressing the free field operator as a sum of plane 
waves weighted by creation and annihilation operators. 
It will be more direct to solve the Klein-Gordon equation 
n-D>(a:) = \I[D^[x) = subject to various conditions. 
One such constraint is the normalization condition pro- 
vided by imposing the canonical value of the equal-time 
commutator [(jj^cj)]. This requires 



l-XD,i.)-D,i.))[__^^-m 



(2.2) 



The zero-temperature solution to the homogeneous dif- 
ferential equation satisfying the above initial condition is 

1 



I i 1 



Att'^ [t + ze)2 - r^ 

At zero temperature D^{x) is analytic throughout the 
lower-half of the complex t plane and D^{x) is analytic 
throughout the upper-half of the complex t plane. For 
later purposes it will be convenient to write these in terms 



^>(^)\t=o'~ 
^<(^)|t=o= 
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The finite-temperature Wightman functions have some 
standard properties [0-0 that are easily demonstrated by 
inserting a complete set of energy eigenstates between the 
operators in (2.1). First, D^{x) is an analytic function 
of complex time in the strip — /3 < Imt < 0, which will 



be referred to as the lower strip. Similarly £)> {x) is an- 
alytic in complex time in the strip < Im t < (3, referred 
to as the upper strip. The values of each Wightman func- 
tions on the boundaries of its region of analyticity are re- 
lated by the Kubo-Martin-Schwinger (KMS) conditions 



11 : 



D<{t + il3,r)^Dy{t,r) 



(2.4) 



The special property of free fields that will be essential 
is that their commutator [0(a;),^(O)] is a c- number at 
all X. The difference of the two Wightman functions is 
the thermal average of this c-numbcr and therefore the 
difference is independent of temperature: 

D>{x) - D^ix) = D>ix)\^^^ - D<{x)\^^^. (2.5) 

This implies that the Wightman functions have the struc- 
ture 

D>{x)^D:,{x)\^^^^ + E{x) 
D<ix)^D<{x)\^^^ + Eix). 

Since Z?> (x) is analytic in the lower strip and Z3< (x) 
analytic in the upper strip, the function E{x) is actually 
analytic in the double- width strip —f3 < Imt < (3. 

Because the fields satisfy the massless Klein-Gordon 
equation, then E{x) does too: \I\E{x) = 0. The most 
general spherically symmetric solution has the form 



E{x)^^(Fi{u)+F2{v)y 



From the definitions in Eq. ( |2.1| ) , time translation invari- 
ance implies that Dy(t,r) = D^{—t,r). Consequently 
E{x) must be an even function of time. This makes the 
two functions arc the same: Fi ~ F2. Therefore the 
Wightman functions have the form 



D>{x) 



D<{x) 
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+ F{v) 
F{v) 



(2.6) 



The KMS condition Eq. (2.4) requires that F satisfy 
1 „, . 1 



+ F(u) = 



F{u~iP). 



(2.7) 



u+ie ' ' u — if3 — ie 

This fixes F(u) to have an infinite number of simple poles: 
1 1 



F{u) = Y. 



u — i{n/3+e) u + i{nP+e) 



(2.8) 



This is the complete answer and has all the correct prop- 
erties. Since Imu — Imi, there are no singularites in the 
closed strip — /3 < Imi < /3 as expected. The nearest 



poles are outside this strip at u = ±i(/3 + e). It is rel- 
atively easy to sum this series and obtain a more useful 
expression. 

Complex time in the open strip: Often one is inter- 
ested in complex time in the open strip 



/3 < Imt < /3, 



which of course includes real time. For the open strip, 
one can set e — >0. Then the series sums to 



F{u) 



ttT coth{TrTu) . 



This has no pole at m = but does have poles at the 
border u = ±i 6 b ecause of the limit taken. Substituting 
this into Eq. ( ^.6[ ) gives 



1 



D>{t,r) = —{S{u)~Siv)) 

--^ — (coth(7rTM) +coth(7rrw)) 

- ^ ( coth(7rTu) + coth(7rrw)) 



(2.10a) 



(2.10b) 



Note that it is only possible for t an d t — z/3 to both be in 
the open strip defined in Eq. (2.9) if both arc complex. 
But if both are complex then the Dirac delta functions 
have no support. Consequently the KMS condition is 
satisfied trivially here. It w ill be satisfied in a nontrivial 
fashion below in Eq. ( 2.14 ). 

Imaginarv time: The imaginary time formalism uses 
t = —IT where r is real. Then u — r — it and v = r -\-iT. 
The two Wightman functions are equal, Z)> (—it, r) = 
£)< {—IT, r) and given by 



D{—iT, r) 



'iT 



sinh(27rTr) 



A-Kr cosh(27rrr) ~ cos(27rTT) ' 



(2.11) 



Naturally this is periodic under t —^ t ± p. 

Arbitrarv complex time: To treat all the poles of F{u) 
in Eq. ( |2.8| ) correctly one cannot use the limit e ^ 0. 
For non-zero e the series can be summed in terms of the 
standard function ipi^z) = dhi[T{z)]/dz to obtain 

F{u) = iTtp[l+T{e-iu)] - iTV'[l+T(e+m)]. 

This is analytic in the closed strip — /? < Imt < /? as 
expected. 

Concise notation: For later purposes it is convenient 
to express F{u) as a derivative. Therefore define a new 
function 



f{u) = T[l+T{e~iu)]T[l+T{e+iu)]. 
Then F{u) = -d\n[f{u)]/dr so that 



(2.12) 



Comparison with Eq. (2.6) shows that the thermal 
Wightman functions can be concisely expressed as 



D>ix)^nd>{x) 
D<{x)=nd<{x). 



(2.9) The lower case functions d{x) are given by 



d>{x) = - 
d<{x) = - 
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fiu)f{v) 



{u — ie){v + ie) 

f{u)f{v) 
{u + ie){v — ie) 



(2.13a) 
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(2.14a) 
(2.14b) 



Note that the difference, rf>(a;) — d^{x), is independent 
of temperature. 

Special case: For many purposes, such as re al ti me, it is 
adequate to use the open strip given in Eq. (2.9), which 
results from e — > 0. In this limit 



lim/(u) = 

e— »o sinh(7rru) 



(2.15) 



This m akes E q. (2.14) rather simple. The ±ie remaining 
in Eq. ( ^.14| ) produce the correct light-cone singularities 
6(u) and 6{v) displayed earlier. In two important cases 
these Dirac delta functions have no support: either t is 
complex or t is real but not on t he lig ht-cone. In either 
case one can omit the ±ze in Eq. ( 2.14 ) in which case the 
two functions are equal: (i> (x) = (i< (x) = d{x) where 



d{x) 



167r2 



■In 



ih(7rTu) sinh(7rTt;)] . (2.16) 



In this regime the two thermal Wightman functions are 
equal: D^{x) — D^{x) =\Z\d{x). 



B. Spin 1 Gauge Bosons 

The same methods can be used to obtain the thermal 
Wightman functions for massless gauge bosons directly 
in coordinate space. In a general covariant gauge the 
Lagrangian density is 

C = -\f,,F^'^ - ^^{d,An\ 

where F^i, = d^Ai, — d^A^ and ^ is an arbitrary gauge 
parameter. The thermal Wightman functions are 

D'i!'{x)^ -iTr{gA^'{x)A''{0)) 
D>iy{x)^ -iTT{gA''{0)A^{x)). 

One can deduce these from first principles using the same 
arguments as employed in Sec IIA. The subsequent dis- 
cussion will show that 

D^''{x) = {-,rn+{l~Od''d-)d>ix) 

D^/ix) = {-v'"'D+{l-Od''d'')d^ix), (2.17) 



where dy{x) and d <;(x) are the functions functions al- 
ready given in Eq. ( 2.14 ). 

(1) The first check of these Wightman functions is that 
they have the correct analyticity in complex time and 
satisfy the correct KMS condition 

D'tr{t-t(3,r)^D>^''{t,r). 

(2) The next check is that the Wightman functions 
must satisfy the correct homogeneous differential equa- 
tion. The variation of the Lagrangian gives 

'^ " - (2.18) 



- P 4'^ 



where the tensor P is linear in the first derivative: 

1 

The equation of motion for the field is 



Po 



Vpfid^ - T]f^^dp 



:VpKdfj,. 



(2.19) 



1, 



= dppp.pA'^ = [ - 77«^n + ii--)d,d^]A^. 

Applying this differential operator to the Wightman 
function in Eq. (^.17) gives 



From Sec IIA, nndy{x) = nDy{x) = and so the 
equation of motion is satisfied. 

(3) The third check is that the field operators in the 
Wightman functions satisfy the correct canonical com- 
mutation relations. The canonical momentum conjugate 
to A'^ is 

n. - Po.^A'^- 

The equal-time canonical commutation relations are 
This requires that the Wightman functions to satisfy 



Pon.iD^'ix) 



D^^^ix)) ^^-S^S'ir). (2.20) 



To check this, use Eq. ( |2.17| ) and ( |2.19| ) to obtain 

Pp.,D^-{x) = (S^dp - S^d, + 7^p.dnD>ix), 

where Dy{x) is the thermal Wightman function for 
scalars. For the appropriate difference of Wightman fmic- 
tions one needs 

p,p,{D^^''{x)~D>:y{x)) 

= iS:do - V^K + VondniDyix) - D^{x)). 

The right hand side contains both time derivatives and 
space derivatives. At t = the spatial derivatives vanish 
because i?>(0, r) = D^{0, r). The above result simplifies 
to 

[Pop,{D'^%x) - D>:y{x))]^^= 5:[{DA^) - D<{x)\^^^. 

Eq. ( ^.2[) determines that the value of the right hand side 
is —ST^ir)^ as required. This completes the proof that 
Eq. ( |2.17 ) is correct. 



C. Spin 2 Gravitons 

Standard quantum gravity is based on the Einstein- 
Hilbert Lagrangian with a gauge-fixing term 

2 ^ 

in which R is the scalar curvature, g = det(gpj,), and 
K^ = 32TrG with G Newton's constant. A conventional 
gauge fixing term is 0,§ 

This corresponds to the Feynman gauge in Yang-Mills 
theories. More general covariant gauges will be discussed 
later. With 77^'' the Minkowski metric, the graviton field 
h^'^ contains all quantum fiuctuations: 



(2.21) 



r^gt^i^ ^ jfi^ + ^ht^i'^ 



Keeping only terms that are quadratic in h produces the 
free Lagrangian density 

Co = l{dph^p){d''h'^^') - \{d,hJ^){dPh/) 
-{dpKp){d''hPP) -f {dphnid^hfso). 
The thermal Wightman functions to be computed are 



A"'"/3/'^~i _ 



fiiyaf3 



(x) = -iTr(£-/t^''(x)/i"^(0)) 
(x) = -iTr(ei/i"'5(0)/i^''(a;)). 



Subsequent argument will show that these are given by 



D 



Mi^a/J/^X _ 



{x)= {~ri'"'T]'"^ + T]''"'q'"^+rit''^T]'''^)D<{x). (2.22) 



That the results are expressed in terms of the scalar 
Wightman functions Dy{x) rather than the potential 
functions dy{x) is a peculiarity of this Feynman-like 
gauge. In more general covariant gauges the Wightman 
functions for gravitons depends on the potentials (i> (x) . 

(1) Because of the properties of the scalar Wigh tman 
functions, the graviton Wightman functions in Eq. ( 2.22| ) 
satisfy 

D'i:'''l^{t-i(3,r)^D'ir''^{t,r), 

and are analytic in the appropriate regions. 

(2) To obtain the equation for the graviton field one 
needs the partial derivatives 



dC 



Pn. 



M'"' 



(2.23) 



where the tensor P is linear in the derivative operator 



-t'^pKX^U \-^^K\pV ^^]K,\^lliy)fJp 



-K. 



• ^ K. ^^K.mi.iJ ^\ 






(2.24) 



and A is given by 



^^hiXpv c) kVk.p'HXu 



VKi^mti)- 



(2.25) 



The differential equation for the graviton field operator 
is d'^ Pp^Xtiuh^'^ [x) = 0, or more concisely 

{-VnxVt^u + Vk^Vxu + Vn^V\f.)nh^"'{x) = 0. (2.26) 

The graviton Wightman functions Eq. ( ^.22 ) automat- 
ically satisfy \Z\Dy^°' {x) — since the scalar functions 
satisfy nDy{x) =0. 



(3) The third check of Eq. ( p.22D is that the gravi- 
ton operators satisfy the correct canonical commutation 
relations. The canonical momentum conjugate to /i'''*' is 

The canonical equal-time commutation relations are 



-z[n.,(a;), h'^'^iO)]^^^ = ~{5:5^ + 5,^5^)5' 



^'r) 



Therefore the Wightman functions should satisfy 



[Po. 



Xp.li 



'D 



pva.f3 



(x)- 



-^piyaf} I 



'D'i!'"^{x))]^ 



iK5^ + K'KW{^ 



(2.27) 



To check that the Wightman function Eq. (2.22) satisfies 
this, first apply the differential operator in Eq. (2.24): 



P. 



pnXfiu 



-^piyaf} I 



',Df'"''{x) 



-j/^I/Q/3^ 



Dt^^^ix)) 



= (S^S,^ + Si5^)dp{D^{x)-D^{x)) 

-{K^x + 5^'5^)d^{Dy{x)-DAx)) 
-{S:6f + S^^S;)dx{D>{x)-D^ix)) 
+iVpn Sf + 77pA<5/)a" {D> (x) - D^ (x)) 
+iVp.S^"+Vpx6:)d^{D^{x)-D^{x)). 

Now set the index p = and t he tim e t = 0. The first line 
obviously coincides with Eq. (2.27) and it is straightfor- 
ward to check that the remaining four lines will always 
vanish. For example, if all of the free indices k, A, /z, v 
are spatial, the last four lines vanish. This, because the 
spatial derivatives on the right hand side are zero since 
Z?> (0, 7^ = Z3< (0, r) . If three of the indices k,X, fi,!^ are 
spatial and one is 0, there will be one non- vanishing time 
derivatives among the last four lines but it will be multi- 
plied by a tensor that vanishes. For the remaining cases 
it is easy to enumerate the possible values of the free in- 
dices K, A, fi, V and verify that the canonical commutation 
relation Eq. (2.27) is fully satisfied. 



General Covariant Gauge : The r esults for the Wight- 
man function given in Eq. ( 2.22 ) are for a particular 
gauge analogous to the Feynman gauge. A more general 
gauge-fixing term is 



^gf- = ^^/M-(3a/=5.9"^)(a/3/=^/''). 

where ^ is arbitrary. The differential equation for the 
graviton field h^^ (x) is then 

{-Vk-XVpi^ + Tli^pVXu + riKuriXii)Uh^'' 

= {'i^-7){Vi^pdxd^ + ri:^i,dxdp + rixpdt,d„+rix^df,dp)h'^'". 

Following the same procedures as above one can show 
that the thermal Wightman function is 



D 



t^^O'P („\ — 



1 



{x)= (- -7;^''77"'3+77^"r;''^+77^'^77'''3)nd>( 



+ {l-j){-2{ri^"'d°'d^^ + r^°'^^^'^'')dy{x) (2.28) 

Note that the general Wightman function d epen ds on all 
partial derivatives of (i> {x) , whereas Eq. ( |2.22| ) , in the 
Feynman- like gauge, depends only on Dy{x) = n'^>(2;)- 

III. FERMIONS 

The methods employed in the previous section carry 
over to fermions except that the KMS condition con- 
tains an additional minus sign. This change results in 
a completely different asymptotic behavior in the deep 
space-like region. 

A. Spin 1/2 Fermions 

The thermal Wightman functions for spin-1/2 fields 
are given by 

S>af3{x) = -iTT{gipa{x)ipp{0)) 

S<ap{x) = iTr(gi'(/'/3(0)V'a(2;))- 

The ordering of the spinor indices is important, but the 
indices will be suppressed in the following. The relative 
sign difference between the two functions is conventional. 
At zero temperature the Wightman functions are 






1 



i 1 

47r2 (i + ieY - ^^ 



The former is analytic in the entire lower-half of the com- 
plex t plane and the upper is analytic in the entire upper- 
half of that plane. 



At nonzero temperature the regions of analyticity are 
reduced. In the complex time plane, S^{x) is analytic in 
the lower strip —(3 < Imi < and 5<(a;) is analytic in 
the upper strip < Imi < /3. On the boundaries of these 
regions the values are related by the KMS conditions M : 



5<(t + */3,r) = -5>(i,f). 



(3.1) 



The difference between the Wightman functions is the 
anticommutator of the fields. For free fields this anti- 
commutator is a c-number: 



Syapix) - S^apix) = -iJV'q (a;) , "0/3(0) } 



(3.2) 



and thus the right hand side is independent of tempera- 
ture. This implies 

S>{x) — 5'<(a;) = S^{x)\rp^f^ — S^{x)\rp^Q. 

so that the Wightman functions have the structure 

S>ix)^Sy{x)\^^^+E{x) 
S<{x)^S<{x)\^^^+E{x). 

The free-field equation —ij-d-ip = implies that E{x) 
must satisfy —iydE = . To solve this, let the unknown 
function E = ij-dF{x) where \ZiF{x) — 0. The general 
spherically symmetric solution is 

E{x)^i-fd-(Giiu) + G2{v) 
r \ 

From their definitions, the Wightman functions satisfy 
[S'>(a;)]^ — 7o<5'>(— a;*)7o and this makes the two func- 
tions the same: Gi = G2- Combining this with the 
previous zero-temperature results allows the full thermal 
Wightman function to be expressed as 



S>(x) 
S<{x) 



{ij-d)ay{x) 
{ij-d)(j<{x), 



in which the new functions have the form 
1 „, . 1 



a>{x) = 
cr>(a;) = 
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u ^ le 
1 



u - 



le 



G{u) 
G{u) 



V + le 
1 



V - 



le 



G{v) 
G{v) 



(3.3) 

(3.4a) 
(3.4b) 



All the temperature dependence is contained in the sin- 
gle unknown function G. It w ill b e determined by the 
fermionic KMS condition Eq. (3.1), which requires that 
ay{t—il3, r) = — cr<(t, r). For the function G this requires 
that 



1 



u — i{(3 - 
The solution for G is 



+ G(u -if3) = - (^— + Gill)) . 
\u + le / 



G(«) = ^(-l) 



1 



1 



u~i{nl3+e) u + i{nl3+e) 



(3.5) 



This is analytic in the closed strip —(3 < Imt < /3. The 
nearest poles are just above this strip at r + t = i{f3 + e) 
and just below the strip at r + t = — i(/3 + e). The alter- 
nating signs will produce more rapid convergence than in 
the bosonic case. 

Complex time in the open strip: For many purposes 
one is interested in either in real time or in complex time 
in the open strip 



- (3 <Imt< p. 

For t in this open region one can set e ^ in Eq. 
which allows the sum to be easily performed: 



(3.6) 



Giu) 



ttT 
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u sinh(7rru) 



M) 



(3.7) 



There is, of course, no pole at u 
u = ±inP because of the limt e - 
and <7<: are 



: 0. There are poles at 
•0. The results for (7> 



(7>{x) 



1 

Sirr 

iT 

Svrr 

1 



{Siu)~Siv)) 

1 1 



sinh(7rru) sinh(7rru) 



^<{x)^—{-S{u) + Siv)) 
STir 

_ iT 
Sttt" 



sinh(7rru) sinh(7rru) 



(3.8a) 



(3.8b) 



It is worth noting that the KMS condition ct> {t — iP, r) = 
~a^{t^ (3) is now satisfied in a trivial manner because of 
the restriction to the open strip in Eq. ( |3.6|) . Because 
of this, it is only possible for t and i — i/3 to lie in the 
open strip if both are complex and if both are complex 
then the Dirac delt a fu nctions have no support. In the 
full solution Eq. (3.5) the KMS condition is satisfied 
nontrivially. 

Imaginary time: The imaginary time formalism uses 
t — —IT where r is real. As occurred for bosons, the two 
Wightman functions are equal, cr>(— zt, r) = (7<(— ir, r) 
and given by 



< 



-iT sinh(27rTr) cos(7rrr) 
27rr cosh(27rTr) - cos(27rTT) ' 



(3.9) 



This is antiperiodic under t ^ t ± j3. 

Arbitrary complex time: Wi thout approximation one 
can sum the series in Eq. (|3.5[) to obtain 



G{u) 



TIT T T 

T ,.1 T , .. T ,^ T . 

Since V'(^) = d\n\r{z)\/dz this is equivalent to 



Giu) = -^Hg{u)]. 



in which the lower-case function g is 



r[i + |:(e + m)]r[i + |(e-zu)] 



r[l + 4-(e + m)]r[l + ^(e-m) 



(3.10a) 



Putting this together in Eq. (3.4) gives the general result 

(3.10b) 



cr<{x) =ns<(a;). 



in which the lower-case functions are 



s>ix) = - 



s<ix) = - 



167r2 



■In 



'liu)9iv) 



167r2 



■In 



{u — ie){v + ie) 

9iu)g{v) 
[u + ie){v — ie) 



(3.10c) 



The next section will show will express the Wightman 
functions for gravitinos in terms of these functions. 



B. Spin 3/2 Gravitinos 

The massless Rarita-Schwinger field Q plays an impor- 
tant role in supergravity as the spin 3/2, supersymmetric 
partner of the graviton. In that context it is referred to 
as the gravitino. As occurred in the bosonic cases, the 
passage from the lower spin fermion (1/2) to the higher 
spin fermion (3/2) is easily accomplished. The free La- 
grangian density for the field i^fi{x) is pO| 

The gauge fixing is necessary so as to break the invariance 
of the first term under transformations ipp —^ tpp + dpTp, 
where -0 is any spin-1/2 field. A convenient choice for 
gauge-fixing is 

The free Lagrangian can be rewritten as 

C = -V^%.V", (3.11) 

where the tensor Lfj_i, is linear in the first derivatives: 

L^"" = ~iij^'d''+Yd^ + V^''l-d-^j''i-fd)Y)- (3.12) 

Each Wightman functions is a 4x 4 matrix in spinor 
space and a rank 2 tensor in the Lorentz indices: 



S:^''{x) = -iTr(gV^(a;)i/'''(0)) 
Snx)^iTT{g^"{0)i^^'{x)). 



(3.13) 



The spinor indices have been suppressed. The relative 
sign difference in the definitions coincides with the spin 



1/2 convention. As in the bosonic case the most effi- 
cient way to proceed is to display the answer immediately 
and then perform checks. In this gauge the Wightman 
functions can be written in terms of the same matrix- 
differential operator L^" 



5^"(x)^L^"a>(2;) 
5r(x) = i^''a<(x), 



(3.14) 



where a:^{x) and cr<(a;) are the basic spin 1/2 functions 
given in Eq. (3.10). The subsequent discussion will con- 
firm these results. 

(1) The spin 1/2 functions cr>{x) and cr<(a;) guaran- 
tee that the gravitino Wightman functions in Eq. ( 3.14 ) 
satisfy the KMS conditions 



5r(t-*/3,r)= 



-^r(i,r-) 



'S'^" 



it, 7^, 



(3.15) 



and have the correct analyticity properties. 

(2) The equation of motion for the gravitino which 
follows directly from Eq. ([3.11 ) is 



Lc^Vix) = 0. 
Therefore the Wightman functions should satisfy 

'{x) = L^^S^'ix) = 0. 
This is easily satisfied because of the identity 

La^L'^'' = -6>:n (3.16) 






and the fact that no'> = and D (t< 



0. 



(3) The third check of Eq. (3.14) is that the gravitino 
field operators must satisfy the canonical commutation 
relations at equal time. For this it is convenient to ex- 
press the partial derivative of the Lagrangian as 



dL 



d{dPtp^) 



ip r^Ap, 



where T,yXp is the matrix 



T^xp = ih,.VXp + IWyp - ^h\lp - 2»Ta)- (3.17) 

Note that Lij\ — —TyXpOP. The canonical momentum 
conjugate to 'ilj'^{x) is 



iix{x)^ij r 



i/AO- 



These momenta obey the canonical equal-time anticom- 
mutation relations 

~i[r{x),nx{Q)]^^^ = -5}:5\r), 
so that the Wightman functions must satisfy 

(5^ (x) - S^^\x))T,xa] ^^^ = -K (53(f). (3.18) 



To check that these are satisfied by Eq. (3.14), first mul- 
tiply by the matrix T to obtain 

For canonical momenta, set p = and the time t = 0. 
Because cr>(0,r) = cr<(0,r) only the time derivatives on 
the right (i.e. n — Q) are nonvanishing. This reduces to 



(5r(a;)-5^''(:E))r,Ao 



t=o 



= 6^[CT>{X) - CT^ix)]^^^ 



Thus Eq. (3. IS) is satisfied because of the properties of 
the functions cr>(x) and (T<(x) from Sec IIIA. 

General Covariant Gauge: The results for the gravitino 
have been displayed in a particular gauge. As was done 
for the graviton in Sec IIC, one can add a more general 
gauge-fixing term than used above. The Wightman func- 
tions can still be expressed in terms of two derivatives of 
the functions cr>(a;) and (J<:{x) in Eq. (B^). 



Deep time-like region: In the deep time-like region de- 
fined b y \t\ — r 3> l/(7rT), the asymptotic behaviors of 
Eqns. (2.1C) and (|3.8|) are 



D>{x) 



-iT 



2Trr 
iT 



"^^1*1 sinh(27rTr) 



o->(a;) -^ ;^e"^^l*l sinh(27rTr) 
ZTrr 



(4.5) 



It is not obvious why both Wightman functions should 
fall exponentially in the time-like region and with the 
same exponent. 

All the Wightman functions for higher spin are simply 
related (by derivatives) to these basic functions. It thus 
appears that in comparison with zero temperature, the 
quantitatively large effects of massless particles at non- 
zero temperature comes predominantly from the T/r be- 
havior of bosons at large spatial separation and from the 



light cone t 



Subsequent publications will explore 



the physical consequences of these results. 



IV. DISCUSSION 
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It is rather surprising that the large distance effects of 
massless particles are, in some sense, simpler at T > 
than in vacuum. At zero temperature the Wightman 
function for spinless bosons is 



D>ix) 



1 



T=0 



47r2 {t - ie)2 - r^ ' 
and for spin 1/2 fermions is Sy{x) — i^^^di^i<jy{x) 

t W * 1 



47r2 [t - ze)2 - r2 



(4.1) 
where 

(4.2) 



In perturbative calculations, the slow fall-off at large dis- 
tances produces long range correlations both in the deep 
space-like and the deep time-like directions. At T > 
the situation is very different as indicated below. 

Deep space-like region: Fo r r — |i| :» l/(7rT) the 
asymptotic behavior of Eqs. ( 2.10| ) is 



D>{x) 



-iT iT 



-27rTr 



47rr 27rr 



cosh(27rri) 



(4.3) 



This behavior can also be understood from the Mat- 
subara formalism, ^-^ which has discrete frequencies 
LOn = 2mTT. The static, n = mode is responsible for 
the leading T/r behavior. Each higher mode produces an 
exp(— a; „7-) fall-off. For fermions the asymptotic behavior 
of Eq. $m is 



ay{x) 



-iT 



2iTr 



-'^^'' cosh(7rTt) 



(4.4) 



Since fermion have odd frequencies, w„ — (2n + l)i:T, 
the exponential fall-off also comes from exp(— [x)„r). 
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